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A concrete proposal for the realization of a Bose-Einstein condensate of kicked rotators is pre- 
sented. Studying their dynamics via the one-dimensional Gross-Pitaevskii equation on a ring we 
point out the existence of a Lax-pair and an infinite countable set of conserved quantities. Under 
equal conditions we make numerical comparisons of the dynamics and their effective irreversibility 
in time, of ensembles of chaotic classical-, and BECs of interaction-free quantum-, and interacting 
quantum kicked rotators. 

PACS numbers: 03.75.Nt, 95.10.fh, 02.30.1k 



The mechanical system known as 'kicked rotator' is 
an important and well studied minimal model describing 
the onset of chaos in Hamiltonian dynamical systems. 
Named after it's original proponents it is also known as 
'Chirikov- Taylor' map and reads |lj 

Pn+i = Pn + Ksinq n+1 

q n+ i = (q n +Pn) mod (27r) (1) 

The map is generated by the Hamiltonian Hr — p 2 /2 + 
K cos q J2 n 5{t — n) by taking stroboscopic cross-sections 
of the dynamics immediately after each kick, p, q is the 
pair of action angle variables of the system. We use units 
in which the moment of inertia 8 and the kicking period 
T are set to 1. Then the only parameter of the system 
is K. For K > 0.976 ... all KAM-curves in phase-space 
are known to be destroyed [2( and the system is chaotic. 
This shows up in a diffusive increase of (p 2 ) with n for 
typical initial ensembles in the (q,p)-plane localized in p. 

This system has also been studied extensively in quan- 
tum mechanics, both theoretically El and experi- 
mentally 0. The quantum system is found to mimick 
the classical diffusive increase of (p 2 ) only for a finite 
time interval |3|, Afterwards the apparent discrete- 
ness of the local quasi-energy spectrum of the quantum 
kicked rotator asserts itself in keeping the expectation 
value (p 2 ) bounded. A discrete local quasi-energy spec- 
trum implies quasi-energy eigenstates, which are local- 
ized in angular-momentum, i.e. square summable over 
the quantized angular-momentum HI £|. This localiza- 
tion of quasi-energy eigenfunctions implies the boundcd- 
ness of (p 2 ) for all n since an initially localized state has 
significant overlap only with quasi-energy states localized 
in its neighborhood. The upshot is that the resulting dy- 
namics of the quantum system is no longer chaotic and 
instead classifies as quasi-periodic. 

An experimental investigation of the quantum kicked 
rotator was achieved via its quantum optical realization 
0. There, mixed-state ensembles of laser-cooled non- 
interacting atoms placed in a periodically pulsed opti- 
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cal lattice could be used to follow in detail the classical 
and quantum dynamics of noninteracting kicked rotators. 
However, strictly speaking the quantum optical realiza- 
tion is only made for kicked particles whose configura- 
tion space R is the covering space of the configuration 
space R/Z = § 1 of the kicked rotator. Classically the dy- 
namics of the kicked particle can simply be projected on 
the configuration space R/Z and is then indistinguishable 
from the dynamics of the kicked rotator. Quantum me- 
chanically the quasi-momentum hk of the non-interacting 
kicked particle with < k < 1 is conserved and equi- 
distributed in the initial ensemble. It can be eliminated 
by a linear shift in the momentum p, however only at the 
cost of modifying the periodic boundary condition of the 
wave function to il>k(<l + 27r) = e l2 ipkil)- 

As a consequence one has realized an ensemble of 
different kicked rotators with boundary conditions la- 
belled by the continuously distributed quasi-momentum 
hk. The total momentum (p + hk) of the kicked par- 
ticle is not discrete in this ensemble, while the quan- 
tized angular momentum p of the kicked rotators is not 
separately accessible to observation. Averaging over the 
quasi-momentum is therefore unavoidable and neccessar- 
ily modifies any observed quantum fluctuations. 

The advent of Bose-Einstein condensates offers new 
possiblities: 

-The new technology to trap and tightly confine ul- 
tracold atoms in optical traps now permits to think 
about new schemes for a direct realization of the quan- 
tum kicked rotator, without the detour via the quantum 
kicked particle. 

- By using an initial state generated by a Bose-Einstein 
condensate it is now possible to prepare an ensemble of 
atoms in a pure state. If the condensate is sufficiently 
diluted in the initial state (e.g. by coherent expansion 
prior to inclusion in the ring-trap) the subsequent dy- 
namics is still that of interaction-free atoms. However, 
distinct from the previously realized uncondensed mixed- 
state ensemble of ultracold atoms, the dynamics of the 
interacting-free BEC will now have to display the full co- 
herent fluctuations of the quasi-periodic dynamics of the 
quantum kicked rotator. 

- Taking advantage of the natural interaction in Bose- 
Einstein condensates the new case of an interacting Bose- 
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running wave, pulsed periodically 

FIG. 1: Proposal for a BEC of kicked rotators. A BEC in 
a 1-d ring-trap of radius R in the (x, j/)-plane is illuminated 
by a periodically pulsed strongly detuned running wave laser 
in z-direction with intensity I engineered to I — Io + -Ti-fj. 
Length of arrows proportional to I. 



Einstein condensate of quantum kicked rotators can also 
be realized. 

Even prior to the modern interest in Bose-Einstein 
condensates this system has been studied as an interest- 
ing model system Q. Recently it has also been studied 
from the point of view of Bose-Einstein condensation Q . 
In particular numerical evidence was given that for suffi- 
ciently strong interaction the quasi-periodic dynamics of 
the quantum-kicked rotator becomes chaotic. 

Here it is our aim to make a proposal for the physi- 
cal realization of the kicked rotator and a corresponding 
Bose-Einstein condensate with optically trapped ultra- 
cold atoms, to derive an infinite number of explicitely 
time-dependent conserved quantities of the associated 
Gross-Pitaevskii equation, and to compare the dynam- 
ics of the three incarnations of the kicked rotator (clas- 
sical, quantum without interaction, and quantum with 
interaction) under otherwise equal conditions. 

For a physical realization (see fig^l we imagine a Bose- 
Einstein condensate enclosed in a quasi one-dimensional 
ring-trap (in the (x, j/)-plane, say). This can be achieved 
e.g. by optical trapping of the condensate in an opti- 
cal doughnut mode. The angle q of the rotator will be 
the coordinate s along the ring, q = s/R, where R is 
the radius. The kicks can be realized [j| by the periodi- 
cally pulsed excitation with a strongly detuned running 
wave laser transverse to the ring (in the z-direction). For 
the cos g-dependence of the kicking potential, which is 
proportional to the laser-intensity in this situation, the 
intensity of the laser beam should be designed such that 
it increases linearly in one direction (say the ^-direction 
with x = R ■ cos q) in the plane of the ring-trap while be- 
ing constant across the trap in the two other directions. 
Thus the periodically pulsed Rabi-frequency fl depends 
only on x, il 2 = fl 2 + flf cosg, and on time t. In a way 
typical for BECs the mean energy must be inferred from 
the expansion of the cloud of atoms after the trap has 
been switched off. 

Now we turn to the dynamics of the kicked rotator 
for the interaction-free classical case, the interaction-free 
quantum case and the interacting quantum case. The 



first two cases have, of course, been extensively studied 
in the literature, and are only given here for the sake of 
comparison. 

The equations numerically solved in the quantum case 
are 
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-V 2 V + Kcos(g) Y, S(t-n)tl>+g\fl>\ 2 il>, (2) 



with g = 0, and g — 5, respectively. Here time is mea- 
sured in kicking periods T, ~% = HT/Q with = mR 2 
is a dimensionless parameter whose size is a measure of 
the importance of quantum effects. The dimensionless 
interaction strength is given by g = AttN Rka/Q where 
Q is the effective cross-section of the trap and a is the 
s-wave scattering length. The kicking strength is given 
by K — (1/A) Jp ulse Qi(t)dt, the time-integral over the 
short pulse in one period, where A is the detuning. It 
should be noted that the spatially constant ilp-part of Q 2 
has no physical effect, since it just adds a constant overall 
phase-shift periodically in time. (In our numerical plots 
we have chosen the three independent dimensionless pa- 
rameters K — 5 well inside the domain of classical chaos, 
% = or 1, and g — or 5. In all cases we choose as 
an initial state the one with p = and q equidistributed 
over the periodicity interval.) 

Before we look at specific numerical solutions we wish 
to point out a rather remarkable general mathematical 
structure of (0 with a surprising additional link to the 
nonintcracting quantum kicked rotator. This structure 
reveals itself by the observation that eq.(0) is the one- 
dimensional non-linear Schrodinger equation in an exter- 
nal time-dependent potential (see e.g. appendix D of 
[Toj')- Setting either the external potential or the inter- 
action to zero this equation is integrable. Remarkably, 
even for nonvanishing interaction and external potential 
the equation still admitts a Lax-pair, i.e. it can be rep- 
resented as the compatibility condition 
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of the pair of equations — Xv,~k^ — T v with 2x2 
matrices 

^ ) V V&Kq,t) ik(q,t) 

T(n f sf A (q,t) B(q,t) 

~\C{q,t) -A(q,t) 
Here k(q, t) can be chosen in the form 



k(q,t) = —ikdhnp(q,t)/dq 



(4) 



where tp(q, t) satisfies the evolution equation of the non- 
interacting quantum kicked rotator, eq.Q with g = 0. 
The matrix-elements A, B, and C of T satisfy three cou- 
pled linear inhomogeneous first order partial differential 
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equations in q, with coefficients, depending on the ex- 
ternal time-dependent potential in (J2J and i/j,ip*,<p and 
their first derivatives, which are all 27r-periodic in q. In 
these equations for A, B, C the time t appears only as a 
parameter. The detailed explicit form of these equations 
will not be needed here beyond the fact that they admit 
solutions A, B and C which are 27r-periodic in q. 

From the Lax-pair, whose existence is often a sign of 
integrability, a countably infinite set of conserved quanti- 
ties can be derived: The compatibility conditions © for 
X and T implies the relation t(d t X ~d q T) + [X, T] = 

Defining by a path-ordered integral the 2x2 matrix 



<E(t> 



U(t) = Pexp 



1 



dqX(q,t) 



(5) 



which is a functional of the condensate wave-function 
xp(q,t), ip*(q,i), and the external explicitcly time- 
dependent kicking-potential at fixed time t, and using 
the 27r-periodicity of T(q, t) in q, T{2n, t) = T(0, t), it is 

possible to show that 



T(0,t),U(t) 



(6) 



This implies the existence of infinitely many (explic- 
itly time-dependent) global constants of the motion 

Cn({i>(<l,t),ip*{q, t)},t) = const with 



C n =TrU n (t). 



(7) 



Unfortunately, we have so far not been able to determine 
the Poisson-brackets of these conservation laws. Conclu- 
sions about the integrability of eq.(0) can therefore not 
been drawn from this analysis at this point. At any rate, 
the dynamics described by (JSJ is highly constrained by all 
these conservation laws. The question whether the dy- 
namics nevertheless is indeed chaotic for strong coupling, 
as numerical results of ; 8J suggest, therefore remains a 
nontrivial one. It is hard to settle numerically since, as 
we will see in the following, the numerical error for the 
strongly interacting system is difficult to controll over any 
sufficiently long time-interval. Of course this fact by it- 
self also lends support to the conjecture that the strongly 
interacting system may be chaotic, the existence of the 
conservation laws notwithstanding. 

Let us now turn to our numerical comparison of the 
three cases outlined before. In the quantum case the 
expectation value of the energy of the kicked rotator im- 
mediately after the n'th kick can be expressed by the 
formula 



K 



E 

m — 1 



dqp m {q)v m {q) sing 



(8) 



This expression is valid for g = and g ^ 0. Here p m (q) 
and v m (q) are the probability distribution and the an- 
gular velocity of the BEC of kicked rotators in config- 
uration space immediately after the m'th kick. They 



FIG. 2: Mean energy {E) of the standard map eq.Q without 
(solid line) and with (dotted line) time reversal for a kick 
strength of K = 5.0. 
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FIG. 3: Mean energy of the kicked rotator, according to the 
GP-equation for g = 0, integrated with the split-operator 
method |l2J. Time reversal (dotted line) t/j — > -0^ is applied 
between the 250th and 251th kick 



are defined by the macroscopic wavefunction ipmiq) = 
y/p m {q) exp(z$ m (g)), v m (q) =-fcV$ m (q). Expression © 
can be directly compared with the corresponding very 
similar formula for a classical ensemble of (noninteract- 
ing) kicked rotators (E n ) = K X^m=i Vm sing m where the 
horizontal bar denotes the average over the initial ensem- 
ble. 

In fig|21 the result for the ensemble of classical kicked 
rotators is shown. This figure is similar to a correspond- 
ing result in [TT| . Along with the well-known diffusive 
increase of the mean energy over 500 kicking periods the 
result of a sudden time-reversal after 250 kicks is shown. 
Due to the chaoticity of the classical map in conjunction 
with the necessarily finite accuracy of the numerical sim- 
ulation the dynamics is effectively irreversible. In figEI 
the corresponding result for the BEC of kicked rotators 
is plotted for the case where the interaction is negligible. 
Dynamical localization leads to the replacement of the 
classical diffusive increase of the mean energy by errat- 
ically and quasi-periodically oscillating interferences of 
the coherent matter wave, which are not averaged out 
due to the purity of the initial wavefunction. Further- 
more, the time-reversal after 250 kicks leads to a com- 
plete retracing of the dynamics to the initial state . 

In fig0]we finally plot the mean energy for a BEC of in- 
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FIG. 4: Same calculation as for the kicked rotator of Fig. [3] 
but with interaction g = 5.0. 

teracting kicked rotators. It is clear that the interaction 
has a strong quantitative but a much lesser qualitative 
influence on the interferences of fig[21 It may be appro- 
priate to stress at this point that figs and 0] were calcu- 
lated with identical numerical codes, the only difference 
being the interaction constants g — for fig|3]and .9 = 5 
for fig^fl (In both cases the angular momentum quantum 
number I was restricted to the interval —64 < I < 63, and 
the 2-7r-interval of q was also discretized into 128 points 
appropriate for the fast Fourier transformation in a 64- 
Bit computation. 16,665 iterations of the split-operator 
method were carried out between two subsequent delta- 
kicks.) It is also apparent from figQ]that the interact- 
ing dynamics becomes effectively irreversible again, even 
though the effective irreversibility can be seen to be much 
weaker than in the classical case. A similar remark ap- 
plies to the diffusive increase of the mean energy. While 
it is clearly absent in fig|3]it reappears in fig01 however 
in a much weaker form and also in a more intermittent 
manner than in the classical case of fig[21 which shows the 
strong influence of the coherence of the nonlinear matter 
wave. An increase in numerical accuracy (which must 
of course remain finite) only increases the time-interval 
over which the dynamics is reversible without changing 
the situation qualitatively. The effective irreversibility 



and the diffusive increase of the mean energy both be- 
come stronger if the interaction strength g is increased. 
This is because the g\ip\ 2 -term introduces a rapidly vary- 
ing g-depcndcnt effective potential, which tends to wash 
out interference effects. An increase of ~fc works in the 
opposite direction, because the more pronounced quan- 
tization of angular momentum will enhance interference 
effects. 

In conclusion, we have considered a Bose-Einstein con- 
densate of kicked rotators. We pointed out the exis- 
tence of infinitely many explicitely time-dependent global 
conserved quantities of the underlying one-dimensional 
Gross-Pitaevskii equation on a ring constraining the dy- 
namics, and gave a comparison of the mean energy of the 
Bose-Einstein condensate with and without interaction 
and the corresponding classical ensemble. For fixed nu- 
merical accuracy we also compared the time-reversibilty 
of the simulation. To achieve a decent numerical accu- 
racy, a large number of time-steps inbetween the kicks is 
required in the case of nonvanishing interaction, where 
the dynamics cannot be reduced to a simple quantum- 
,map (see also 8] ) . As we showed here by the reversibil- 
ity check this severely limits the number of kicks over 
which the dynamics is quantitatively reliable and hence 
reversible. It is tempting to interpret the lack of re- 
versibility as another sign of chaos as in the classical case, 
but the existence of the Lax-pair and the conservation 
laws seem to call for caution. The experimental realiza- 
tion of the BEC of kicked rotators, possibly along the 
line suggested here, both in the non-diluted interacting 
and in the diluted non-interacting case would be very in- 
teresting and might shed more light on the nature of the 
dynamics. 
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